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Abstract 


Second-order effects upon the wave systems associated with 
bodies in uniform supersonic flight are considered. An appropri- 
ate scaling of the independent variables is introduced to take 
into account the locally small, but cumulative, nonlinearities 
which undermine the validity of the classical perturbation theories 
at large distances from the body. 

For planar flows, the solution valid at large distances from 
the body (the wave structure solution) is uniformly valid to the 
body surface, and the boundary conditions may be applied directly. 

For this case, the simple wave results of Friedrichs as corrected 
and extended by Lighthill are confirmed, with a minor correction 
noted for the location of the rear shock at very great distances 
from the body. It is noted that effects locally of third order 
must be taken into account behind the trailing shocks to give proper 
global integrals relating to the lift and drag of the airfoil sec- 
tion and to determine the positions of the trailing shocks to second 
order. 

For flows about finite bodies, the wave structure solution must 
be matched with a quasi-cylindrical , local solution to provide the 
inner boundary condition. This procedure gives rise to an inter- 
mediate (or one-and-one-half-order) solution which corresponds to 
the first-order solution over a slightly modified body. In these 
two wave structure solutions, the azimuthal angle enters only as a 
parameter, and the flow in each azimuthal plane may be thought of 
as that about some equivalent body of revolution. Dependence upon 
azimuthal angle does arise in the true second-order wave structure 
solution. The complete general solution has been found for this 
second-order equation, but the third-order local solution is required 
to effect the matching which determines the inner boundary condition. 
For flows over bodies whose surfaces lie everywhere near an axis 
aligned to the flow, an additional matching is required with a solu- 
tion valid near the body surface to determine the local solution. 

This s lender-body solution is also required to third order to de- 
termine the wave structure solution to second order. The lack of 
a local second-order theory of sufficient generality is the major 
hurdle in determining the wave structure solution. Local third-order 
solutions can presently be obtained for only the simplest conical 
geometries . 

Predicted shock angles for the flow over slender wedges (planar 
flow) and cones (axisymmetric flow) are compared with exact inviscid 
calculations to show the improvement afforded by inclusion of second- 
order effects. Inclusion of the full second-order solution is re- 
quired to achieve appreciable improvement over the first-order theory 
for the conical case. 
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Chapter One 


The Introduction 


l.A. General Introduction 


The present treatise is concerned with the mathematical descrip- 
tion of the wave systems which emanate from bodies in supersonic 
flight, and, in particular, the behavior of these systems at large 
distances from the body. A perturbation theory, which includes ac- 
count of essential cumulative nonlinearities, is considered, accurate 
to terms of second order. 

Such a problem is of practical interest as it relates to the cal- 
culation of sonic boom overpressures, where we are interested in pre- 
dicting the effects of a high flying, supersonic aircraft in the vi- 
cinity of the ground. It is also of theoretical interest as it illus- 
trates techniques useful in solving problems of weakly nonlinear wave 
propagation. For, while the classical linearized theories of super- 
sonic aerodynamics are perfectly adequate for the prediction of dis- 
turbances near the body (and especially surface pressures), the failure 
of these theories at large distances from the body has long been recog- 
nized, and a correct theory must in some way account for cumulative 
nonlinearity within the wave system. 

The equations of steady, compressible fluid flow are inherently 
nonlinear, even when the simplifications of negligible viscosity and 
irrotationality are used. The mathematical treatment of nonlinear 
equations is still largely restricted to especially simple cases, and 
the complexity of the fluid flow equations precludes treatment in any 
generality without additional assumptions to provide further simpli- 
fication. 

Perhaps the most useful simplification -- at least for the case 
of supersonic aerodynamics -- is to consider that the flow is every- 
where a small perturbation upon some known flow field (usually uni- 
form) . The nonlinear problem is then broken up into a hierarchy of 
simpler (often linear) problems, and the solution may be found as an 
asymptotic series, successive terms of which contain factors of higher 
order in the perturbation parameter. 

The classical approach to perturbation theory for steady super- 
sonic aerodynamics yields the scalar wave equation for the velocity 
potential, to be solved subject to appropriate boundary conditions, as 
a first approximation. Higher approximations also satisfy the wave 
equation, but with additional nonhomogeneous terms which are known 
functions of the lower-order solutions. Thus, a hierarchy of linear 
problems is set up and may be solved to any degree of accuracy desired. 
Since the homogeneous part of each equation is always the scalar wave 
equation, the characteristics are not revised in successive approxi- 
mations. Higher approximations may then be thought of as giving suc- 
cessive terms in the Taylor series expansion of the exact solution 
about the linearized (freestream) characteristics (Van Dyke (1952)). 


Such a scheme has been extremely useful for predicting distur- 
bances near bodies immersed in supersonic streams, especially for 
predicting surface pressures. However, since the characteristics 
are not revised in higher-order solutions, the predictions become 
increasingly inaccurate as the distance from the body is increased. 
Such behavior is a result of the "nonuniform validity" of the lin- 
earized solution. I.e., the solution obtained by such a scheme is 
not equally accurate at all distances from the body. Physically, 
this is due to the fact that the actual (exact) characteristics and 
the linearized (freestream) characteristics are getting farther apart 
as they are followed away from the body. This deficiency is serious 
if we are concerned with accurately describing the flow field at 
large distances from the body, and is interesting from a theoretical 
viewpoint as it illustrates the nature of (at least one type of) non- 
linearity which must be accounted for to achieve uniformly valid 
solutions . 

The essential nonlinearities in the solution at large distances 
from the body may be accounted for in two ways. We could solve the 
purely linear problem for the solution on the linearized character- 
istics, and simultaneously correct the placement of these character- 
istics in physical space due to the perturbations of the solution. 

Or, we may try to rescale the independent variables of the problem 
in such a way as to bring the essential (or cumulative) nonlineari- 
ties into the lower-order problem. Hopefully, this rescaling will 
also provide some simplification in the other terms of the equation, 
so that the new equation will be mathematically tractable. In this 
second approach, the first-order equation will be nonlinear (as it 
must be to achieve our goal), but higher-order equations will again 
be linear. 

This second method will be used here. The specific purpose of 
the present treatise is to discuss the nature of the second approxi- 
mation, or the second term in the asymptotic series representation 
of the perturbation solution, valid uniformly to large distances 
from the body. 

The purpose of such a second-order theory is two-fold. First, 
it may provide increased numerical accuracy for solutions simple 
enough to be calculated analytically. But, secondly, and more im- 
portantly, a study of the second-order solution is one way to verify 
the uniform validity of the first-order solution, thus showing that 
the entire procedure is indeed a rational approximation to the exact 
solution in some limit. Further, such a study gives greater insight 
into the nature of the first approximation, and provides limits on 
the accuracy of that theory. 

l.B. Relevant previous researches 

The background of the linear theory of supersonic aerodynamics 
is well established in the literature and textbooks, and will not be 
discussed here. (See, e.g., Ward (1955) or Heaslet and Lomax (1954)). 

Local second-order effects were considered by Van Dyke (1952) . 

He solved the potential equation by iteration (i.e., using the first- 
order solution to evaluate quadratic terms) to obtain a second 
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approximation. He found complete particular integrals for the non- 
homogeneous terms in the cases of planar and axisymmetric flows, thus 
reducing the second-order problem to an equivalent first-order prob- 
lem. His solutions gave marked improvement for the calculation of 
surface pressures, but were not uniformly valid at large distances 
from the body. He recognized this fact, and realized that it was due 
to the inability of the theory to revise the characteristics in suc- 
cessive approximations. 

Friedrichs (1948) first obtained solutions accurate to second 
order and uniformly valid to large distances for the planar flow over 
% airfoils, using the fact that the flow is (to second order) a simple 
wave. Correction of his expression for the shock slope in terms of 
the perturbation quantities on either side and inclusion of effects 
due to the broad third- order waves behind the trailing shocks (which 
have a cumulative effect of second order) were discussed by Lighthill 
(1954). 

Whitham (1952) rendered the first-order solution for the flow 
over slender, axisymmetric bodies uniformly valid to large distances 
by applying the coordinate straining technique of Lighthill. This 
essentially uses the linear solution to correct the placement of the 
characteristics in the physical plane. 

Hayes (1954) used the method of rescaling the independent vari- 
ables to order terms according to their cumulative (rather than local) 
effects to obtain the proper form of the first-order wave structure 
for more general flows. For flows about finite planar systems, the 
method relies upon the method of "coincident signals" of Hayes (1947) 
to supply the appropriate boundary condition in terms of the large- 
distance asymptotic behavior of the local solution. 

Work currently in progress by Landahl et al, (1968), and some- 
what parallel to the present study, has looked briefly at second-order 
corrections to the wave structure theory of Hayes. They have inde- 
pendently noted the need for an intermediate homogeneous wave struc- 
ture solution for flows about finite systems (which we term of one- 
and-one-half order) , which is determined by the local second-order 
solution. This result is used to support the need for experimental 
surveys in the absence of a complete second-order local theory. 

The nonuniform validity of the linear theory at large distances 
is equivalent to the singular behavior of the linear solution near 
the Mach cone in problems of the conical type. This is because the 
relevant parameter is the ratio of the distance a characteristic has 
travelled from the body to the axial distance along the body from 
which the characteristic started, and this ratio becomes very large 
for waves near the leading Mach cone. This nonuniformity may be rec- 
tified using either of the two methods previously mentioned for the 
wave structure problem. Lighthill (1949) applied the method of 
strained coordinates to the problem to determine the shock wave po- 
sition to first order (i.e., its lowest-order deviation from the 
freestream Mach cone). Broderick (1949) used a separate expansion 
scheme valid near the Mach cone (essentially the method of matched 
asymptotic expansions) to prove that his expression for the second- 
order surface pressure on a circular cone, derived without consid- 
eration of the shock, was correct to that order, and also verified 
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Lighthill's expression for the position of the shock to first order. 
Bulakh (1961), in the Soviet Union, showed Broderick's result for a 
general conical flow, and again verified Lighthill's result for the 
shock wave position to first order. 

The current work also bears a formal resemblence to the second- 
order transonic theory cf Hayes (1966). The effects which exert a 
cumulative influence upon the supersonic wave system are locally im- 
portant in transonic flows. Thus, the similarity coefficients of 
the two theories (in the sense of Van Dyke (1958) as giving the form 
of the solution) are essentially the same. 

l.C. General nature and scope of present theory 

The present treatise presents a perturbation theory for steady, 
inviscid supersonic flows which is uniformly valid to large distances 
from the body of interest. The theory is basically an extension to 
include second-order effects of the first-order wave structure theory 
of Hayes (1954b) . Uniform validity of the solution to large distances 
is obtained by rescaling the independent variables of the problem such 
that nonlinear terms are ordered according to their cumulative (rather 
than local) effects in an outer, or wave structure, region. The solu- 
tion so obtained is uniformly valid to the body surface for planar 
flows, but the inner boundary condition for flows about finite systems 
must be obtained by matching with an appropriate local solution. 

The complete theory of planar flows is presented, along with a 
sample solution. The results obtained confirm those of Friedrichs 
(1948), obtained by a characteristics method, as corrected and ampli- 
fied by Lighthill (1954), except for the position of the rear shock at 
great distances from the airfoil, where the Lighthill result is in 
error (see Appendix B) . 

For the case of flows about finite systems, the matching with the 
local second-order solution reveals the necessity of a one-and-one-half- 
order wave structure solution, satisfying homogeneous equations. This 
corresponds to a first-order solution over a slightly modified body. 

The third-order local solution is required to obtain the second-order 
wave structure solution. The complete general solution of the second- 
order wave structure problem has been found, but difficulty in obtain- 
ing the third-order local solution to provide the inner boundary con- 
dition precludes calculation of any examples (except the rather spe- 
cialized case of a right circular cone at zero angle of attack) . 

In the following chapter, the basic equations are developed, and 
a brief discussion of the scaling of the wave structure region is in- 
cluded. Chapter Three specializes the analysis to planar flows. The 
general form of the solution is given, as well as some specific fea- 
tures noted from the sample solution. Flows about finite systems are 
treated in Chapter Four, including the matching problem with the local 
solution. The most important results are summarized in Chapter Five, 
along with some general . remarks on future researches. Details not 
essential to the development of the text are included in the appendices. 
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Chapter Two 


The Development of Basic Equations 


2. A. Potential Equation 

As the heading of this section suggests, we shall work almost 
exclusively with a velocity potential. This allows us to describe 
the vector velocity field by a single scalar equation, but not with- 
out the loss of some generality. That such a potential exists de- 
pends heavily upon our assumption of small perturbations, and par- 
ticularly that there be no strong shocks in the flow. For, while 
the flow is initially irrotational , and considered to be inviscid, 
we cannot completely disregard the entropy increments caused by 
shock waves, and it is vorticity arising from the variation in 
strength of these discontinuities that would eventually force us 
to abandon the use of a velocity potential if we were interested 
in a theory more accurate than the second order. 

It can be shown that, for our purposes, a velocity potential 
exists everywhere in the flow field to the required order except 
behind the trailing shocks, where entropy and velocity (pressure) 
perturbations are of the same order. (See Appendix A). This is 
essentially because the entropy increment at a weak shock wave is 
of third order in the strength of the wave, and we are concerned 
with only second-order effects. We are interested in the third- 
order quantities behind the trailing shocks only because they have 
a cumulative influence upon the position of these shocks which is 
of second order (See section 3.C.). 

We consider a steady, homocompositional , homoenergetic, in- 
viscid fluid flow. The equations of continuity, momentum conser- 
vation, and isentropy for such a flow may be written 


v-(^) = ° , 

+ -i-v-p = O , 

vS= o , 


( 2 . 01 ) 

( 2 . 02 ) 

(2.03) 


where ^ is the fluid density, is the vector velocity, ^ is 
the pressure, and S' is the specific entropy. 

We form a modified continuity equation by converting the den- 
sity gradient to a pressure gradient, using the condition of isen- 
tropy : 



(2.04) 
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where <X_ is the acoustic velocity. 

Then, eliminating the pressure gradient between this and the 
momentum equation (2.02) yeilds 



Or, introducing the velocity potential. 


$ 


(2.05) 


such that 



(2.06) 


we have 


= V§- y('a'7§ v ) . 


(2.07) 


This is the fundamental equation for all steady, irrotational 
flows. It is an exact equation, subject only to the requirements 
of steadiness, homocompositionality , homoenergy, and homoentropy. 

It is, however, nonlinear and extremely difficult to treat with 
any generality. (Note that the nonlinearity arises not only from 
the right-hand-side of the above equation, but also from the fact 
that Gt* is also a function of the potential, $ ). 

Therefore, additional assumptions must be made to cast the 
problem as one that is mathematically tractable. A standard tech- 
nique (and the one we sha}l use) is to assume the flow may be des- 
cribed as everywhere a small perturbation upon a uniform streaming 
motion. This allows us not only to expand the right-hand-side of 
the equation, neglecting higher-order products of the perturbations, 
but also to expand the acoustic velocity in powers of the pertur- 
bations. Thus 

a* = f fS ) } ( 2 . 08 ) 

where P and are thermodynamic parameters of the fluid (always 
understood to be evaluated in the undisturbed stream) and defined 
by 

r - ±(^X , 

^ = “^(fps • 


(2.09) 

( 2 . 10 ) 
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The symbol £ designates the difference of the quantity it precedes 
from its value in the undisturbed stream, and the subscript ( 3^ 
refers to evaluation in the undisturbed stream. For a calorically 

perfect gas, P - .1^! and <r = O . 

In order to use this expansion in the potential equation, we 
must have a relation between the pressure perturbations and the 
velocity (potential) perturbations. Such a Bernoulli relation may 
be derived as follows. For a homocompositional, inviscid flow, the 
Crocco-Vaszonyi theorem states 


= VH - Tv$ } (2.11) 

where H is the total enthalpy and T is the temperature of the 
fluid. This may be combined with the momentum equation to give, for 
a homoenergetic (i.e., H = constant) flow, 

which may be expanded and integrated to give 

Jf - ^F-f- r^fS * Offs: SptS -) , 

where U is the undisturbed, or free-stream, velocity, 
also write a formal, thermodynamic expansion of this integral. 


( 2 . 12 ) 

(2.13) 
We may 


And comparison of the last two results, by an iterative process, 
yields 



(2.15) 


This is the required Bernoulli equation which relates pressure 
perturbations to changes in velocity potential. Here we have 
introduced the perturbation potential, <$) . 


(j{ £* + , 


(2.16) 
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where *r is a measure of the maximum body slope, c is the body 
length (chord) f and with . We have also in ? 

troduced the Prandtl-Glauert coordinates 

W = Vc , 

> (2.17) 

? - /* Z /c , 

in place of a cartesian system (x,y,z), with its x-axis aligned to 
the freestream, and have used the notation 



The Bernoulli equation may 'be substituted back into equation 
(2.08) to give 

of- (2.i 8 ) 

Finally, this may be substituted back into the fundamental poten- 
tial equation (2.07), the right-hand-side of that equation expanded 
in terms of the perturbation potential, and we arrive at 

* 2(r-0 (2.19) 

+ (l-<r) <f)f <}>ff 

-t v<f)\ l <fi)} +‘ m \ 

This equation for the perturbation potential, 0 , is the 

starting point for all perturbation theories, and in the above 
form is accurate to third order in the local perturbation quan- 
tities. It is an extension of the equation given by Van Dyke 
(1952) to include all third-order quantities (not just those im- 
portant near a slender body), and is derived for a general fluid 
(not a calorically perfect gas). 

2.B. Shock relations 


A peculiar feature of nonlinear partial differential equations 
is the presence of regions in the physical plane where the solution 
is multiply defined, due to the overlap of characteristics coming 
from different points on the boundary or initial data. In our case, 


8 



this difficulty is overcome by the insertion of shocks * or discon- 
tinuities in the solution (or its derivatives), across which the 
solution jumps from one branch to another in the characteristic 
plane. The intermediate portion of the solution is then omitted 
to render the solution unique in the physical plane. 

These shocks are inserted according to the classical Hugoniot 
conditions for conservation of mass, momentum, and energy fluxes 
across the discontinuity surface, and the requirement that the 
tangential velocity be continuous across the surface. 

We describe the shock location as 


« • > • (2 ' 2o) 

and use square brackets, [ ] , to denote the jump in any quantity 
across the discontinuity surface. Thus, the requirements of con- 
tinuous tangential velocity and normal mass flux may be written 

= ° > ( 2 . 21 ) 

and 

( 2 . 22 ) 

respectively, where Q* is the velocity component along the normal 
to the discontinuity surface. The conservation of energy is auto- 
matically assured through our use of the Bernoulli equation (2.15) 
to relate pressure, velocity, and entropy perturbations. Finally, 
in place of the usual momentum relation, we shall use a quantitative 
expression for the entropy jump 

■pTS] = + O(C^), (2.23) 

which is valid for weak shock waves. (See, e.g., Hayes (1954a.).) 

We may expand the density as 



where . 


yU * i (4r-o-- 7 ) 


(2.25) 
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and 


l) 



(2.26) 


and these last two parameters are always assumed to be evaluated 
in the freestream. For a calorically perfect gas, j^ = f- 

and "0 * T - f . 

Substituting from the Bernoulli result (2.15) gives 


f tWj- * C r-o»‘6f} 

gXc -^ r - ?) } 

-!^>T-5S *•••), 


(2.27) 


and since 


f' 


•/ 1 + /s^r s f ' 


(2.28) 


we may write (2.22) as 

[- T/3(0 f + W>-V^} - riV{-t('flp L - <j> f 

+ ^’{(i-dhV, **) - 

+( 0 f C fc]' 0 . 

The term involving the entropy jump may be eliminated using (2.23) 
to give 


P-V 

V/->VE£ 


• (0 f+ vf^) + ) J 

-(^orwy * O. 


(2.30) 
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This last equation and that expressing continuity of the potential 
across the discontinuity, (2.21), will be used to determine the 
locations of any shocks in the flow field, and the jumps in per- 
turbation quantities across them. 

2. C. Boundary and Initial Conditions 

The appropriate boundary condition for inviscid flows is that 
the velocity normal to the surface at every point on a physical 
body vanish. In addition, we require the initial condition that 
all perturbations vanish upstream of the body. (In linearized 
theory, this initial condition is applied on the leading Mach cone; 
we shall apply it immediately upstream of the leading shock sur- 
faces, which are slightly ahead of the leading Mach cone.) These 
conditions form a well posed problem for the wave-type equations 
with which we shall be dealing. 

For planar flows, the wave structure solution is uniformly 
valid to the body surface, and the boundary condition may be applied 
directly. The mean surface approximation is used, and the conditions 
are applied along the line = O . 

For flows about finite bodies, the boundary condition is used to 
determine a local solution, which is valid in the vicinity of the 
body. The large-distance asymptotic behavior of this solution is then 
matched with the small-distance asymptotic behavior of the wave struc- 
ture solution to determine the potential in this outer region. Match- 
ing of the potential and not merely its radial derivative (which might 
be though to be a sufficient condition) is required to completely de- 
termine the positions of any shocks in the flow. 

2.D. General nature of scaling concepts 

The appropriate scaling of variables in the wave structure region, 
i.e., that region in which cumulative nonlinearities must be taken in- 
to account, may be seen in several ways. For simplicity of exposition, 
we will illustrate the arguments only for the case of planar flows. 
Similar arguments may be followed for more general flows, but the es- 
sential ideas are more easily lost in algebraic complexity. A discus- 
sion of the scaling problem for supersonic airfoil theory is also in- 
cluded in Van Dyke's book (1965), to which the following paragraphs 
owe. 

The most physical argument is one based upon the hyperbolic nature 
of the equation describing the flow, and the failure of successive 
approximations to correct the characteristics to more closely approxi- 
mate those of the full nonlinear equation. Since the straightforward 
perturbation procedure attempts to correct the solution by expanding 
it in a Taylor series about the linearized characteristics, the pro- 
cedure will become inaccurate when the distance between the actual 
characteristic and the linearized characteristic starting from the 
same point on the body becomes of order unity, for then the first 
correction is of the same order as the linear approximation. Since 
for the plane flow over a slender obstacle of thickness of 0(-r) the 
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perturbation in signal velocity (or characteristic velocity) is of 
0(-r) and constant along the characteristics in the linear solution, 
the slope of the exact characteristic is everywhere of 0(-t) differ- 
ent from that of the linearized one passing through the same point. 
Thus, if we follow the linearized characteristic a distance of 0(1/t) 
away ' the body, we will arrive at a point in the flow that is of 

0(1) . from the point we would have found by following the exact 

charav Stic. This suggests that if we rescale the lateral vari- 
al 1 -- e. , the one which measures distance along the characteristic 

-- . xch lat it is of 0(1) when the physical variable is of 0(l/-r), 
the equation should rearrange itself to provide an accurate descrip- 
tion of this long range, cumulative nonlinearity. 

A second argument (which, of course, arrives at the same result) 
is more formal in nature, and is seen by observing the nature of the 
nonuniformity in the straightforward perturbation solution. For the 
planar flow over a bump described by 


j b = TC-P(Vc) \ 


(2.31) 


the straightforward expansion gives for the perturbation potential 


0 = 4> 0 + + ■■■ , 

the solutions (Van Dyke (1952)): 

= - -PC * -pfi > 

f *“7 5 

# - - + f- of {f'to} ^ 


f' 1 - P% 


(2.32) 


(2.33) 


(2.34) 


Now, the solution (2.32) is a uniformly valid asymptotic expansion 


as r-*> 0 if ^ <f> { 


grows no faster than 0 O in any region. How- 
ever, due to the presence of the third term in the solution (2.34) 
for <p , the contribution from the second-order term will be of the 


same order as that from the first-order term when 


fcr 


This , 


V " r 

too, suggests that in the region of 0(l/-r) distant from the body, 
a new scaling of variables is required. 
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Chapter Three 


T he Theory of Planar Flows 


3 . Equations and Boundary Conditions 

For planar flows, we assume the flow properties to be indepen- 
dent of ^ , and introduce the contracted, semi- characteristic 

coordinates 


where 



(3.01) 


(3.02) 


is the fundamental lateral scale factor which forms the basis for 
Hayes’ (1954) first-order similitude. The new variables, F and 
V| , correspond essentially to a phase and an age variable, re- 
spectively, in a geometric analysis of the analogous one-dimensional 
unsteady problem. (Or, they may be thought of as the short- and 
long-time variables in a multiple time scale analysis.) 

The equation for the perturbation potential (2.19) then becomes 

V = f{j\ + 2 ( + V 

+ 4 (r- 1) $ T <t> t («--<) <b*<f>^ 

We introduce the expansion for the potential 


} 


(3.03) 


1 > * 4 



(3.04) 


whence we arrive at the hierarchy of equations 


’V'M’rt = °' C3 ' 05) 

VV (3 - 06) 

+ 4 (r-/)& f <£ 0 ^ + (cr-04ol<t>o ff . 
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We note that the first-order equation (for <f> 0 ) is nonlinear. 

In addition to the first term (which is essentially the wave equa- 
tion in semi- characteristic coordinates for waves of primarily one 
family) , there is now a term to allow correction of the character- 
istics. This latter term is, of course, absent in the linear theory. 
The second-order equation is again linear (as would be the case for 
all higher approximations) with nonhomogeneous terms which are known 
functions of the first-order solution. The extreme simplicity of 
the first-order equation is a major factor in our ability to obtain 
useful results. 

We expand the shock shape in a manner similar to that for the 
potential, 


= H 0 (v|) + ; 

the shock relations (2.30), (2.21) become 

(3.07) 

[*,+*; - ■ ° > 

(3.08) 

[$*,+ = o ; 

(3.09) 


and 


l<t>, + 

+ J- Q ^o r f= 0 > 


(3.10) 


C \ h,^ + + F. 7 * rf )3 - ° i t 3 - 11 ) 

where (2.21) has been differentiated along the shock surface to 
produce (3.09) and (3.11), and all quantities are evaluated at 
the surface E© , in analogy with the mean surface approximation 
used for the boundary conditions. The differentiation of (2.21) 
has the effect of introducing an arbitrary constant into the ex- 
pression for the shock location. This constant is eliminated in 
the first-order result by starting the shock at the first point 
of characteristic overlap. In the second-order result, (3.11) 
is used only to simplify (3.10), and the nondifferentiated form 
of (2.21) (eqn. (3.18)) is used to locate the shock. 

The first-order shock relations (3.08) and (3.09) may be 
combined to give 
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and 


_ -L j _ ~T~ 

a L %-] ^°5 > 

L^ + J = ° • 


(3.12) 


(3.13) 


The first result states the well known fact that (to first order) 
the shock velocity is the mean of the characteristic velocities 
on either side of the discontinuity, and the latter may be inter- 
preted as stating that the Riemann invariant of right-running 
(returning) waves is unchanged in crossing the shock. (See below). 
This latter result, combined with the first integral of the first- 
order equation (3.05), 


> (3.i4) 

and the fact that all perturbations must vanish everywhere up- 
stream, gives the relation 


+ Z - ° , 


(3.15) 


which will prove most useful in simplifying the second-order 
problem. 

The second-order equation (3.06), with the use of (3.15) and 
(3.05) may be simplified to 


4 r? +C4 f * f ), = 


(3. 16] 


The second-order shock relations (3.10) and (3.11) may be 
combined to give 

CVW ■ R + 

ry _ CM*- . i)/ ~\r* t3 C3.17) 

where the result 


a * j -o , 


(3.18) 


15 



obtainable directly from (2.21), and the identity 


* §L^ 3 3 - (3.19) 

have been used, as well as the simplifications of (3.05) and (3.15). 
The identity (3.19) may be derived from (3.09) by multiplying within 
the brackets by <p 0 , and using the fact that 

Ofe>] = a. } (3.20) 

where the bar denotes taking of the mean of the quantity on either 
side of the discontinuity. Equation (3.17) essentially gives the 
change in the Riemann invariant for the crossing waves, and (3.18) 
will be used to determine E, 

The wave structure equations presented here correspond to an 
orderly expansion of the characteristic relations for steady, in- 
viscid flow. These equations state that 

-^- d r ±Jd =o ( 3.2 i) 

along 

41" = 1^(6*^^) > (3.22) 

respectively . In the above form, the relations are valid for a 
rotational flow. When expressed in terms of our wave structure 
variables, the characteristic equations become 

4 - j * O (3.23) 

along 

^ * V + I } + ' • • , C3. 24) 

and 

(3 25) 

-j - 0 

along 
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(3.26) 



Thus, the shock relations (3.13) and (3.17) may be interpreted 
as giving the change in the Riemann invariant of the crossing sys- 
tem of characteristics. This change is zero to first order, but is 
non- zero in the second approximation, due to the inability of the 
shock to match perfectly with the conditions in a simple wave. This 
locally third-order effect, when followed to large distances from 
the body, contributes a cumulative effect of the second order. From 
(3.17) and (3.25), the change in the invariant, Q, for the crossing 
wave is 

HQ!! - 2&P 1 - * ^12^) J . (3.27) 

This quantity will play an important role in the determination of 
the location of the rear shock at very large distances from the 
airfoil. (See Appendix B) . 

Since the flow under consideration is two-dimensional, the wave 
system is nearly planar not only in the wave structure region, but 
near the body surface as well. This fact is manifested mathematic- 
ally in the uniform validity of the wave structure solution in the 
region of the body surface, with the result that the boundary con- 
ditions may be applied directly to the wave structure solution, 
without the necessity of introducing an inner (or local) solution. 

We consider the shape of the body surface to be described by 

TC-f(f) ) (3.28) 

where f(f) is the body shape function, and the subscript ( )^ refers 
to the body surface. Then the appropriate boundary condition is 
that the streamline slope, 

-foe = -T +•■ , (3.29) 

be tangent to the surface at every point on the body. 

Since we are considering bodies which lie everywhere near the 
mean plane vi = 0 > this condition may be satisfied at the surface 
by an appropriate Taylor expansion about values on that plane. 

(This is the so-called "me an- surface approximation.") Thus, sub- 
stitution of (3.04) gives 
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(3.30) 


= -f'('E) , 

% “ * ■f'V ‘ > C3 - 31) 

as the required boundary conditions, to be applied at ^ = O . 

To summarize, we here collect the equations, shock relations, 
and boundary conditions in their final form for the first-order 
problem. 



0 ^ 

-—3 

4* 

ii 

o 


(3.05) 


+ = o > 


(3.13) 




(3.12) 


- -f (?) ( on 

7 * ° ; 

(3.30) 

for 

the second-order problem. 



S 


(3.16) 

°V 

f 7 = f" ^ * 6/5* ~ ^ * 


(3.17) 

r, 

-- W V W ,3 , 


(3.18) 



/>*,*• 

■ 2 yS’-r<» f J o* y »0 . 

(3.31) 
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3.B. General solution 


The first-order equation (3.05) has characteristics described 

by . 



along which 0 Of is constant. Thus, the characteristics are straight 
lines in the 5 -^ plane. The value of on a particular char- 

acteristic is determined from the value of f' at the value of g 
corresponding to We denote this value of E as . 

Thus the solution may be written 



(3.33) 


where 

= ° (3<34) 

determines implicitly. This is equivalent to Whitham f s result, 

which states that the linear solution is correct even at large dis- 
tances from the body if the linearized characteristics are replaced 
by the exact ones (or at least by ones accurate to first order). 

The above result reduces to the linear result near the body (i.e., 
for small ), where approaches E and the linearized and 
first-order characteristics are indistinguishable. 

If we are interested in the actual value of the potential, we 

have 

$>= -fCZi.) Zi *1 f(\) > (3.35) 


where the relation (3.13) has been used to eliminate the arbitrary 
function of ^ arising from the integration. The locations of 
any shocks which may appear in the flow are then determined by in- 
tegrating (3.12) in any region where the characteristics overlap 
in the physical plane, starting at the point where any such overlap 
first occurs. 

The second-order equation has the same characteristics as the 
first-order equation, and along these the homogeneous part of the 
second-order potential is constant. Integrating the second-order 
equation (3.16) once with respect to E yields 

+ s (f , (3.36) 
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where S(^) is an, as yet, arbitrary function, introduced by the in- 
tegration. Since the velocity perturbations must vanish everywhere 
ahead of the leading shock, S(^) must be identically zero in that 
region. Indeed, since it is a function only of ^ , it would have 

to be zero everywhere were it not for the fact that it can be in- 
cremented at shocks. This is seen from the shock relation (3.17), 
which states 

W 7 - V'r J C3,i7) 

Thus, while the function $0|) cannot depend upon E in any inviscid 
region, it must be allowed to increase by 



at each shock wave in the flow field. The quantity 


«L + 5^6, 

is related to the Riemann invariant of the wave family crossing the 
shock, and the amount by which it changes is a measure of the de- 
parture of the wave system under consideration from a simple one. 
Thus, the equations (3.36) and (3.17), by comparison require 




(3.37) 


where the summation is over all shocks upstream of the point of 
consideration at any value of 

Now, the first-order solution, <f> 0j , is a function only of Tf* , 
whence 



4 ) T (%) ff fc>J , (3 . 


38) 


and 


where a function of integration to be determined by the 

boundary condition. 
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Differentiating the solution with respect to E , and evaluating 
it at v] = O, we have 



(3.40) 


whence the boundary condition (3.31) requires 





(3.41) 


and we have, finally. 




This gives the complete second-order potential in terms of the body 
shape function, -mT*) 9 in the first-order characteristic coordinate 
system (■£<.> vj ) which is known from the first-order solution. 

The second-order correction to the shapes of shocks in the flow 
is then obtained by evaluating 
curve, aitd using 



e, 



(3.18) 


The uniform validity of the solution may now be seen from a 
comparison of the general solutions for the first- and second- 
order problems. We thus have for large distances that 



(3.43) 


while 



(3.44) 


Thus the ratio of second-order perturbation quantities to those of 
first order' along any characteristic is 



M V-£- rM \ 

\ 2 . + /V t fa) 


(3.45) 
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which approaches zero uniformly as 'T-^O. 


3.C. Sample solution 


The general form of the wave structure solution for planar 
flows is presented in the previous section, and the complete solu- 
tion may be constructed once the body shape function, -ffct,), is 
known. This has been done in Appendix B for the particularly 
simple example of a symmetric, circular-arc airfoil at zero angle 
of attack, and some particular features of the solution will now 
be presented. 

Since the flow will be symmetric about the chord-line of the 
section, only the flow in the upper half-plane is considered. 

Thus, we consider the shape function to be given by 



= O 


> 


-4 < , 

-t < E k . 


(3.46) 


The 
limit of 


first-order solution, 
large “*7 ) is then 


far from the 


airfoil 


(i.e. , 


in the 



(3.47) 


in the region between the shocks described by 





(3.48) 


and is zero elsewhere (i.e., upstream and downstream of the lead- 
ing and trailing shocks, respectively). 

The second-order solution between the shocks behaves, for 
large >7 , as 


+ + , (3.49) 

while the shocks are displaced by an amount 

S, L = ■£( pl + . (3 - 50) 
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The dominant second-order effect at large distances is a dis- 
placement of the shock wave which strengthens the leading wave and 
weakens the trailing wave. These results agree with those of 
Friedrichs (1948) when our analysis is specialized to the case of 
a calorically perfect gas, and his expression for the shock slope 
is corrected as noted by Lighthill (1954). 

Behind the trailing shock, the analysis gives as the poten- 
tial 



which is a function of ^ alone, whence the pressure perturbations 
are of third order in this region. But the form of (3.51) suggests 
that there is a right-running wave that is locally of third order, 
but has a cumulative effect of the second order, and hence must be 
included in the analysis. This effect was first noted by Lighthill 
(1954), and must be accounted for if certain global integrals re- 
lating to lift and drag are to be correct to second order (See 
Appendix C) . This wave eventually interacts with the trailing shock 
below the airfoil (or is M ref lected 11 and interacts with the trailing 
shock above the airfoil in our symmetric case), and alters its po- 
sition to second order. This interaction takes place only very far 
from the body, and we introduce for this "Lighthill region" the new 
scaling 



K z *( • 

Since entropy and pressure perturbations are of the same order behind 
the shocks near the airfoil, the strength of this wave must be deter- 
mined using the roational method of characteristics. In the region 
of interaction with the trailing shock (3.52), however, the entropy 
perturbations are again negligible, and the potential behind the 
trailing shock is 

$*(^*"^+3^//+? f (3.53) 

and the second-order shock displacement is 

— / — rw*- i - jl . rtf- ,1 

X s IjJTzifK " 6p<-- , (3.54) 

which at very, very great distances from the airfoil (i.e., for 
large vf ) , approaches 
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(3.55) 


r 'v 
**•1 


Vff- 


~3Af- + 




This result differs slightly from that of Lighthiil (1954) for 
the position of the trailing shock, which is in error. (See Appen- 
dix B) . 


5.D. Structure of Shock Origin 


Since the characteristics are not revised in our second-order 
approximation, we cannot hope for the theory to give detailed infor- 
mation about the wave structure near the origin of a shock within 
the flow field. The theory does approximate such local behavior in 
a global sense, however, and gives a second-order displacement of 
the point of formation which corresponds to the axial displacement 
of the continuation of the shock shape valid away from the point of 
formation, to the same lateral position as the first-order shock was 
first formed. This is illustrated in Appendix D, where the geometry 
of the true characteristics is investigated in the neighborhood of 
such a point. 

A shock is formed in the flow field when characteristics coming 
from different points on the boundary data first cross. The shock 
then forms a barrier which the characteristics are not allowed to 
cross, thus rendering the solution unique in regions where the over- 
lapping characteristics would otherwise have caused the solution to 
be multivalued. The case examined in Appendix D is that of the shock 
formed in the planar flow over a step formed by first a compressive, 
then an expansive, arc of a parabola. The geometry and first-order 
characteristics are shown in Figure 3.01. The first-order theory 
predicts a shock to be formed with finite strength at (0 } Ajl), up- 
stream of which the flow is undisturbed, and downstream of which 
the solution from points on the body for which 1 is valid. The 

second-order theory then predicts the shock to originate at the 
slightly displaced point, 




(3.56) 


To interpret this in terms of local structure, we must look 
at the geometry of the actual characteristics in the region near 
the point of formation. When second-order effects are included in 
the characteristic relations, the characteristics coming from the 
forward part of the body O < < ] form an envelope described 

parametrically by 


qjT+ 2/0 , 
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Figure 3.01: Geometry of Shock Origin Problem. 


When the coefficient, 

* = £ + 7^- , (3.58) 

is positive, the envelope appears as in Figure 3.02, and a shock is 
formed to intercept the characteristics before the envelope is formed. 
The shock begins with zero strength and the flow is undisturbed ahead 
of the shock. Beyond the point corresponding to , the shock is 

determined by the appropriate relations between the free-stream char- 
acteristics and those coming from the expansive part of the body 
C I ), as in the regular second-order theory. If this slope is 

continued back to , the shift in virtual origin is exactly 

E., = -3X = > (3. 59) 
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whence, the second-order theory is correct in this global sense. 



When the coefficient, X , is negative, the envelope appears 
as in Figure 3.03. The shock begins with zero strength and its 
early development occurs with disturbed flow on both sides. Since 
the shock is now formed for ^ , the virtual origin at V) - ^ 

is merely the value of at that point. It is ' 

K, w = -k* - S, > (3.60) 


again showing the second-order theory to be globally correct. 

When the coefficint, X , is zero, the shock again, as in 
the first-order theory, starts at the point ( Jfc.) , and the 
predicted second-order displacement is zero. 
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Chapter Four 


The Theory for Finite Systems 


4. A. Introduction 


The study of planar flows, while interesting pedagogically , is 
of little practical value, since any body we are likely to conceive 
for supersonic flight will most certainly not look two-dimensional 
from the large distances with which we are primarily concerned. Thus, 
it is important that we should consider the problem of wave structure 
for flows about bodies which are of finite extent (and, perhaps, look 
somewhat like real airplanes). 

The easiest problems of this type to analyze are those of axially 
symmetric flows, since these may still be described with only two in- 
dependent variables. The consideration of axially symmetric flows is 
important for another reason, however, and this is because they hold 
the key to description of the wave structure for flows about more gen- 
eral bodies. This is suggested by the result of Hayes (1947) that, 
for finite planar systems, the perturbations in any azimuthal plane 
look, to first order, like those from some equivalent body of revo- 
lution, if the observer is sufficiently distant from the system. 

(This is the so-called supersonic M area rule 11 , first used for cal- 
culating wave drags). A similar result may be obtained for bodies 
whose every lateral dimension is small (slender bodies), with the 
additional simplification that the Equivalent" body of revolution 
is the same in all azimuthal planes. Thus, since a general body may 
be thought to be made up of suitably arranged slender and planar bodies, 
the result holds for quite general finite bodies. 

Now these results are based upon the large-distance asymptotic 
behavior of the local (straightforward perturbation) solution -- and 
this is exactly what is used to provide the inner boundary condition 
on the wave structure solution. Thus, the first-order wave structure 
depends upon the azimuthal angle only as a parameter, which enters by 
matching with the local solution in the appropriate plane (Hayes (1954)). 
We shall refer to this property as quasi-axisymmetry . 

This property suggests that the appropriate coordinate system in 
which to study the wave systems for finite bodies is a cylindrical one, 
aligned with the freestream. Thus, the equations for the wave structure 
are expressed in a cylindrical coordinate system, with the lateral (or 
radial) coordinate properly contracted to allow consideration of cumu- 
lative nonlinearities. The inner boundary condition is then provided 
by asymptotic matching with the local solution. 

We are hampered in this sense by a lack of local second-order so- 
lutions of any generality. Only for flows over bodies of revolution 
does a reasonably complete second-order theory exist. Thus, we shall 
consider in detail the matching for these flows, and then briefly dis- 
cuss the nature of problems encountered for flows about non-axisymmetric 
and planar bodies. 
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Finally, it will be seen that the local second-order solution 
determines only a one-and-one-hal£-order solution in the wave struc 
ture region. The third-order local solution is required to deter- 
mine the second-order wave structure solution. Third-order local 
solutions can be obtained only for the simplest conical geometries 
at the present time. 

4 . B . E quations and Shock Relations 

The quasi-axisymmetry of the first-order wave structure sug- 
gests the coordinate system 


^ - t : ~ jx . } + S’* > 

^ = 2 . Kt/k , 


( 4 . 01 ) 


for description of the wave systems emanating from finite systems. 
It is essentially a cylindrical coordinate system, aligned with 
the freestream, with the radial coordinate contracted to allow 
consideration of cumulative nonlinear effects. We introduce the 
reduced potential 




( 4 . 02 ) 


which automatically accounts for the geometric attenuation of cylin 
drical waves, and is of order unity throughout the wave structure 
region. 

Thus, the potential -equation becomes 

~ ^77 + ^ + ^7 Vjf) 

+ ‘Kr-0%Vfr + c 4 -°« 

We introduce the expansion for the potential 

*■■■ , ( 4 . 04 ) 

(where the term is included to satisfy the matching require- 

ment with the local second-order solution) and arrive at the set 
of equations 
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'kff + ° > 

+ ' °* 

\ 7 + OW, “ f f^S? * 


(4.05) 


(4.06) 


) 


77 ' " f WJ 7 7 "" (4.07) 

- t( v - V «f> 5 ^ r ct +4 V/ 

ffr 


-i&W' . 


Similarly, the shock relations (2.30) and (2.21), when trans- 
formed as above, with the introduction of 


Vy)= *.(y)+ + E^ 1 e, cy) + ••■ 

^ (4.08) 

become 

L^ + ^4 - ] - 0, 

(4.09) 

[>..+ E^Vv,] - 0; 

(4.10) 

and 

^ V \ + z Vv VV S' 1 "* 

/ ' ^ -1 

' (4.11) 


V ^ \vj + ? 0^\ f 


(4.12) 


and 


30 



- S*^ V*SS,) - jC^\„ 4 Vs, 

- V ^ V 2 \%- 2 N\„ > - &-< (n% t , 

v ^*h. '^’i^’m'^’n* ( 2 - v, i u °’i 

■ v* - ^Ss-^4- ^ t ., x,-m 

* T 1 v ^ v T^V V|\^ - O , 


(4.13) 


C-^*, -- ^,\-< '(.-i ■* ?^.Vi t + 
+ *,0^ + 

* ±r^. nl - + 


(4.14) 


where (2.21) has been differentiated along the shock surface, and 
all quantities have been evaluated on the first-order shock surface, 
tc , by a mean-surface approximation. 

The first-order relations, (4.09) and (4.10) may be combined 
to give 


CV + 


■Kl = c? 


(4.15) 


and 




tv.p 

LK1 


= V- 


(4.16) 


where is the mean of Vby on either side of the surface. This 

latter relationship is recognized as the bisector law, while the 
former may be combined with the first integral of the equation of 
motion (4.05) to give 


¥1^ -+ i: Vo; = O 


(4.17) 
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everywhere in the flow field, in analogy with the planar result. 
This relation will again be extremely useful in simplifying the 
second- order problem. 

The one-and-one-half-order shock relations (4.113 and (4.12) 
may be combined to give 

D^-+ ] = O > (4.18) 

where the result 


- 0 ) (4.19) 


obtainable directly from (2.21) and the fact that h. satisfies 
(4.06) have been used. The former result, (4.18), combined with 
the first integral of (4.06) give the one-and-one-half-order 
analogue of (4.17): 


= 


O 


(4.20) 


everywhere in the flow field. 

Combining the two second-order shock relations yeilds 


0,. + (- 1 2 

A)lS2'- r 


where the result 


(4.21) 


O, + * t <\ 5 )] = o 


(4.22) 


has been used, as well as (4.05), (4.06), (4.17) and the result 

] = - i C^I 3 + -§ LV 1 . (4.23) 

The result (4.17) may be used to simplify the second-order 
equation to the form 
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(4.24) 


-■&OivVk n ) 

■*■ "p y* Ok + ^ # 9») “ i^jA T ^ f 5 • 

For convenience, we here summarize the equations and shock 
relations in final form for the first-order problem: 



(4.05) 

0. ? + *%*;)- o. 

(4.15) 

■^0%| "" Vo^. J 

(4.16) 

-and-one-half-order problem: 

* CWr - ° . 

(4.06) 

t\+W]*°, 

(4.18) 

Zh.-'^VlSI i 

(4.25) 


and for the second-order problem: 

r/y*- 


V ( Wr (- 7+ 7»-*/Wv 
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> 




(4.21) 
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(4.26) 


It is o£ interest to note that as the wave propagates away 
from the body, geometric effects and wave interaction effects re- 
main of comparable magnitude in the second-order equation. This 
is seen from the fact that, at large distances, the wave region 
will be approximately parabolic in the plane; i.e., will 
scale as the square of ^ . (This is because the characteristics 
are straight lines in the plane, and at large distances the 

body appears as a point. Thus the shocks, which define the wave 
extent, are approximately the bisectors of a centered fan and 
parallel lines (the freestream characteristics) ; they are there- 
fore nearly parabolic.) Thus the term 






which arises from the wave interaction, must remain of the same 
order of magnitude as the term arising from the changing geometry, 


Thus, at even very large distances from the body, second-order 
wave structure effects are never dominated by the geometry, and 
both effects must be considered. 

4.C. General Solution 


The first-order solution is almost completely analogous to 
that in the planar case. The equation (4.05) has characteristics 
in any azimuthal plane, O , given by 

^ = \ . (4 - 27) 

along which is constant, whence they are straight lines in 

the azimuthal plane. (Due to the fact that > how- 

ever, they will no longer be straight in the physical plane, but, 
rather, parabolic.) Thus we will again have 

% = - r '( , (4.28) 
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where 


» O , C4.29) 

and Ffrw a ) is a function to be determined by matching with 
the first-order local solution in the appropriate azimuthal plane. 

The potential is then given by 

^ (*<., 7 ; ©) = - 6) + -ky r &) . C 4 .3°) 

Shocks are inserted exactly as in the planar case, by integrating 
(4.16) in regions of overlap of characteristics in the plane. 

The azimuthal angle, © , thus enters the first-order solu- 

tion only as a parameter, and the wave structure in any azimuthal 
plane is the same as that for some equivalent body of revolution. 

This equivalent body is defined by matching with the local solu- 
tion in the azimuthal plane of interest. 

The one-and-one-half-order equation has the same characteristics 
(4.27), and the potential, , is constant along these. The 

equation may be written 



(4.31) 


where the arbitrary function of integration has been set to zero 
by (4.18) and the application of the initial condition. Thus, 


= Q (£t»; ©) > 


(4.32) 


where QO^e) is to be determined by matching with the local 
second-order solution. The correction to the first-order shock 
shape is given by (4.25). 

Again, as in the first-order solution, the azimuthal angle 
enters only as a parameter, and determines the plane in which the 
matching is to be done to determine G 

The characteristics of the second-order equation are again 
given by (4.27), and the homogeneous part of the solution is con- 
stant along them. Integration of the equation (4.24) with respect 
to 3E gives. 


\* %9f vl *- 

♦JjSfyjCv- '■ si ”* 






(4.33) 
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where the shock relation C4.21) determines 





(4.34) 


and the summation is again over all upstream shocks at any value 
of . Thus, the first integral may be written 


rr' 

!*• 


©)= F’fae) +(j£-z)^p. 

_i p ( y1^ F+ 4 ^ Z[o + + <tei] C4 ■ 

+ (" * + §j5*“ > 
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which may be integrated to give 


^ e) = (-f - sjsL rfo; e)^ + (*- j ^ l ) ^ 

+ 4 + ^ I + 4 ^ '^6©^ ) 

-%^- r , + - h & *> . 

The function H (?„; 0) is to be determined by matching with the 
third-order local solution. 

The second-order correction to the shock shape is given by 


(4.36) 


" DJ'Oj] ^ + ^xf) ^ • (4.26) 

The azimuthal angle enters the second-order solution paramet- 
rically in the determination of H(T^6), and also through the 
dependence of F and upon © in the particular integral 

terms of the solution . 

The uniform validity of the finite-body solution at large 
distances from the body may be seen from a comparison of the orders 
of magnitude of the various terms as we follow any characteristic 
to infinity. Thus, for large distances, the first-order solution 
gives 



(4.37) 
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while the one-and-one-half- order solution gives 


L gfo; 

» 1 FTn.e) ’ 


and the second-order solution gives 



(4.38) 


(4.39) 


Thus, it is seen that, in the limit as T-*-0, the one-and-one- 
half-order and second-order contributions to the solution are always 
small compared to that of the first-order solution. Unlike the 
planar case, the higher-order solutions actually become small com- 
pared to the first-order solution at large distances (by a ratio of 
approximately the inverse square root of the distance from the body) . 
Note also that the second-order solution grows logarithmically com- 
pared to the one-and-one-half-order solution. Thus, in a rigorous 
sense, the true second-order solution must be included whenever the 
one-and-one-half-order solution is, although in practical cases, 
logarithmic terms may usually be considered to be of order unity. 

The wave structure solution is not, however, uniformly valid 
for small . This may be seen from the equation for the second- 

order potential (4.24) which has nonhomogeneous terms in reciprocal 
powers of. ^ which will become large for small values of ^ . Or, 

more directly, the nonuniformity may be seen by comparing the second- 
order solution (4.36) with the first-order solution (4.30). This 
comparison shows that the second-order contribution to the potential 
will be of the same order as that of the first-order solution when 


y = 0(t) t 

or, equivalently, 

JE- 0 ( 1 ). 


Thus, in the region 7t = 0(1), a separate expansion must be carried 
out, and the solution of this problem matched with that of the wave 
structure problem. The solution in this region is, of course, the 
classical, "straightforward perturbation" solution. 

The reverse of the above procedure, i.e., an examination of 
the nonuniformity of the local (straightforward perturbation solution), 
is exactly what led to our choice of the contraction of radial co- 
ordinate to be used in the wave structure region scaling (as was done 
in Section 2.D. for the planar flow case). 
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4.D. Matching with Local Solution 


In order to see how the functions rfae) , £(*<.;$) , and 
H (*nG) are determined, it is necessary to look at the local so- 
lution in some detail. The simplest problem which illustrates the 
nature of the matching involved is that of the flow over an axi- 
symmetric disturbance on an infinite cylinder, aligned to the free- 
stream. Thus, if the body surface is described by 


A b = A 0 + > 


(4.40) 


the perturbation potential, expressed as 

4 - & * + " , 


is determined to be 


?-* 


- 






and 




(4.41) 


(4.42) 


(4.43) 


The source functions , and Mm are determined by the first- and 

second-order, boundary conditions 


/ 




and 


(4.44) 


(4.45) 


(See Van Dyke (1952)). The local solution is completely determined 
once these Volterra integral equations are solved f or c?G0 and MCx) . 
This solution may then be used to determine the wave structure solu- 
tion by the matching principle of Lagerstrom which (following Van 
Dyke (1964)) states that 


the m-term inner expansion of (the n-term outer 
expansion) = the n-term outer expansion of 
(the m-term inner expansion). 


We will here carry the matching only to terms of second order 
in the local solution (one-and-one-half-order in the wave structure 
solution), since the third-order local solution is not known. Thus, 
the one-and-one-half-order outer expansion of the second-order inner 
solution is 



Tivfe' 


) 


(4.46) 


where a subscript in parentheses indicates the complete solution 
for the potential to the order indicated, and the subscript on the 
parentheses indicates the region in which the solution is valid. 
Thus, e.g., the subscript (1)^ signifies the solution in the inner 
(or local) region, up to and including terms of second order in that 
region. 

The one-and-one-half-order wave structure solution, expanded 
to terms of second order in the local region is then 


( )./ ~ {- r (*> - KVT ^tr) - t? G (?) ■* - j . 


(4.47) 


Expressing (4.47) in terms of the wave structure variables for com- 
parison with (4.46) gives 





2 . & Mr p 

^ /3 ^ 


(4.48) 


whence the matching requires 



o 


(4.49) 
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(4.50) 



Note, in particular, that the local second-order solution de- 
termines only the one-and-one-half-order solution in the wave struc- 
ture region. To determine the second-order wave structure, we would 
require the local solution to third order, the homogeneous part of 
which would contribute a term like 


2-K 

1 . P* J 


(4.51) 


in (4.46), which would match with a term like 


21T mV* 

i r 


H(rj 


(4.52) 


in (4.48) to determine the homogeneous term in the second- order 
solution (4.36) . 

4.E. 'General Finite Bodies 


The treatment of flows about more general finite bodies is 
severely limited by the lack of a sufficiently general second-order 
local theory. This is primarily a result of the inability to find 
a complete particular integral for the nonhomogeneous terms of the 
equation when non-axisymmetric terms are included. However, a few 
general remarks may be made about these flows. 

a. Quasi -cylindrical flows 

For the flow over a nearly axisymmetric perturbation on an in- 
finite cylinder aligned with the flow, it is useful to think of the 
body surface as described by a Fourier series. To illustrate, we 
will consider only the first term (e.g., the dipole) of such a series, 
though all terms are of equal importance in determining the wave 
structure, since higher-order moments do not decay as in the slender 
body case (see following section) . Thus the body is assumed to be 
described by 

3 % = + TyQj r f(t) + . (4.53) 

The first-order solution will then be of the form 


40 



(4.54) 




r-A 
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- cos© f / f-^ 4 -V(f 1 ^ t - 7 » aL 1 7C \ c? z QQJa 

J 0 V /l ?- A + ^(f-x ) 1 -^ 1 


where ^,(>) and o&(A) are related to the shape functions £ fo) and 
ftOfw’S by Volterra integral equations through the boundary condition, 
which may be applied at Jl 0 The second-order local solution will 

have a similar form for the homogeneous terms, plus particular inte- 
grals of the nonhomogeneous terms which we are unable to calculate: 



O 


ft 


jdjXvbj 


(4.55) 


+ 


( pafticwla** 




ten^O . 


Matching this with the wave structure solution results in 


and 


F (*w,e) 



^ CM ~t 2 . CeS© c?^(a) 



(4.56) 


Q(r b -ei = . 


(4.57) 


The second-order source functions, and >WiWcannot be deter- 

mined, however, since the complete particular integral (to be 
used in applying the boundary conditions) is not known. The second- 
order outer behavior of the first-order local solution matches appro- 
priately with the inner behavior of the wave structure solution, but 
the check cannot be made for the nonhomogeneous terms of the second- 
order local solution. 

As in the axially symmetric case, the local third-order solution 
is required to determine the wave structure solution to second order. 


b . Slender body flows 

For flows over bodies whose every cross-stream dimension is 
small, the following redefinitions are useful. We define the slender - 
body perturbation potential 
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and the new parameter 



rM^r 1 ' 
/ 3 * ' 


whence, with the expansion 


i> = a + ■*••• 


(4.59) 


(4.60) 


the equations of Section 4.B. again describe the wave structure 
when (f is everywhere replaced by W and is replaced by 

In this case, an additional scaling region must be consid- 
ered within the local region, in the immediate vicinity of the 
body. In this slender-body region, radial velocity perturbations 
are large compared with axial ones. The solution in this region 
must then be matched with the local solution, which in turn pro- 
vides the inner boundary condition for the wave structure solution. 
The first-order slender-body solution is a harmonic function in 
the cross plane, plus an additive function of the axial variable. 
Thus, non- axi symmetric terms decay rapidly, and the first-order 
local solution is axisymmetric. The slender-body dipole distri- 
bution contributes to the second-order local solution. Since 
the first-order local solution is axisymmetric, none of the diffi- 
culties of the preceding (quasi-cylindrical) section appear, and 
the wave structure solution, in principle, may be determined (to 
one-and-one-half order). The slender-body quadripole distribution, 
as well as particular integral terms arising from the non-axisym- 
metry of the local second-order solution will enter into the third- 
order local solution, which is needed to determine the second-order 
wave structure solution. 

c. Finite planar systems 

For flows over finite planar bodies, the local theory exists 
only to first order (Hayes (1947)). The so-called n area rule" is 
a large-distance local theory -- which is exactly what we need for 
matching with the wave structure solution. It is based upon the 
fact that the disturbance at a point far from a finite planar body 
is primarily due to sources near the intersection of the forward 
Mach cone from the point of interest and the planar body, and that 
this region may be shrunk to a point, as viewed from the distant 
point of interest. 


For the thickness, or symmetrical, solution, we need consider 
only volume sources, and the first-order disturbance potential at 
a point , is of the form 

0 = -[ [—*: ~ * *** (4.61) 

where the integration is carried out over the region upstream of 
the intersection of the forward Mach cone from the point of inter- 
est and the plane of the body (assumed to be O ) , and is 

proportional to the local downwash due to the thickness distribution. 
The one-term wave structure expansion of this integral may be easily 
shown to be 



1-i 


( 4 . 62 ) 


where 


V = X, - *2. C<75 0 


(4.63) 


and 

o^(V) *»)</**. . < 4 - 64 ) 

Thus, we note that at large distances, the perturbation potential 
from a planar distribution of sources is identical to that of an 
equivalent axisymmetric distribution, given by ofy(v) . Note that 
will depend upon the azimuthal angle, © , through the 

definition of V 

Similarly, for the lifting problem, the potential due to a 
distribution of elementary horseshoe vortices may be represented 
as 


, fCrCs-x,)^ (*i,0 ^ 

yf^ ' <4 ' 65) 

where is proportional to the local lift coefficient. In the 

wave structure region this takes the form 

(<L , •><*». %(»?■!*' . (4.66) 
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In this case, the azimuthal angle enters not only in describing the 
cuts across which the integration will be carried out, but explicitly 
in the integrand. 

Again, the second-order local problem cannot be solved due to 
the lack of a complete particular integral, but it may be noted that 
the homogeneous potential will also be of the form (4.61) or (4.65) 
and will, in the matching with the wave structure solution determine 
the one-and-one-half-order solution Qfrbi©) • 

4.F. Li gh thill region 


The broad third-order wave found in the planar case behind the 
trailing shocks also exists in the finite body case, with a cumulative 
effect of second order. The wave in this case is even broader in the 
physical space, and the Lighthill region defined by 

sr = 

vf - it, , 

corresponds to 

A = 0(V /0 ) 


(4.67) 


(4.68) 


in the physical plane. The strength of the wave cannot be deter- 
mined as it was in the planar case, however, due to the fact that 
the geometry now enters the characteristic equations, precluding 
direct integration. The strength might be determined for specific 
examples, but the great distance at which the interaction with the 
rear shock wave occurs justifies, in some sense, its neglect. 
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Chapter Five 


Summary and Concluding Remarks 


A second-order perturbation theory for the inviscid, supersonic 
flow of a general fluid past a -body has been presented, with partic- 
ular emphasis on the behavior of the wave systems at large distances 
from the body. 

The theory justifies confidence in the first-order theory (which 
is a standard tool in all sonic boom work), by showing that it is the 
first term in an asymptotic approximation scheme, valid in the limit 
of vanishing body thickness, and gives additional insight into the 
nature of that approximation. The theory also provides increased nu- 
merical accuracy for flows about bodies of finite thickness in the 
examples calculated, provided the disturbances are not too large. 

The theory of planar flows generally confirms the results of 
Friedrichs (1948), as corrected and extended by Lighthill (1954). An 
error is pointed out in the latter work, wherein the position of the 
rear shock at very great distances is incorrectly given. It is also 
demonstrated that for shocks originating within the flow field, the 
predicted second-order shift in point of origin corresponds to the 
shift in virtual origin necessary to account for the local structure 
near that point (which structure the present theory does not predict). 

The theory of flows about finite systems is essentially new. Of 
particular importance is the fact that the local solution (i.e., that 
valid ne&r the body, where the cylindrical ' nature of the waves is 
essential) must be known to third order to determine the second-order 
wave structure solution. An intermediate, homogeneous solution (termed 
of one-and-one-half order) is determined from matching with the local 
second-order solution. Numerical results from the calculation of the 
shock position in a conical flow indicate that appreciable increase in 
numerical accuracy is obtained only when the full second-order solution 
is calculated. The one-and-one-half-order solution offers little im- 
provement over the first-order theory. 

The wave structure theory, as presented, is essentially complete 
for the large distance behavior of these systems. Extension of the 
theory to third order would probably not be worthwhile. Such a theory 
would not necessarily give increased numerical accuracy even if its 
algebraic hurdles could be overcome, and the essential nature of the 
higher approximations is illustrated in the second-order theory. 

The most useful avenue for further study is consideration of the 
second approximation for non-axisymmetric bodies in the local region. 
Only the outer asymptotic behavior of this solution is required to 
determine the corresponding wave structure. Determination of the local 
third-order solution is probably not feasible. 
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Appendix A 


The Existence of a Velocity Potential 


An irrotational vector field, i.e., one having zero curl, may- 
be described as everywhere the gradient of a single scalar function, 
called the potential. The irrotationality of the field is a neces- 
sary and sufficient condition for the existence of such a potential. 
We will now investigate the conditions under which the velocity field 
in our particular problem is irrotational. 

Since we assume the flow to be inviscid and initially uniform 
(and hence irrotational) , vorticity may be generated within the flow 
field only at curved shocks, where entropy gradients are incurred. 

The vorticity is related to these entropy gradients for a homocom- 
positional, inviscid flow by the Crocco-Vaszonyi theorem 


|>T » -V" v-S + vH 


(A. 01) 


where 


V 



(A. 02) 


is the vorticity. Thus, for our homoenergetic flow, the magnitude 
of the vorticity may be estimated by 


V ~ 




U 


> 


(A. 03) 


which, using our quantitative expression for the entropy jump at 
weak shock waves, (2.23), becomes 


■ tA ' 04 

Using equation (2.28) for the normal component of the velocity, 
this becomes 


Xt 2 ( 7 


V 


W). 


(A. 05) 


This gives for the magnitude of the vorticity in the planar 

case, 
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(A. 06) 


T ~ 



and for the finite system case. 


V ~ 


cv^jy. 

o 


(A. 07) 


We can now calculate the magnitude of the vorticity in a more 
general Rotational) flow. This is a flow for which a potential 
does not exist, hence the magnitude of the vorticity is given by 
terms of the sort 


These may be estimated by calculating the cross derivative terms in 
the potential analysis. The vorticity is then found to be of order 

K 3 rl7 

— (A. 08) 


in the planar case, and of order 


(A. 09) 


in the finite body case only in the third-order solutions. 

Thus, the velocity potential may be used with no loss in gen- 
erality for describing these flows to second order. 

From the above, it is seen that vorticity must be taken into 
account in the third-order problem, and a velocity potential can- 
not be used to describe the third-order wave behind the body. How- 
ever, in the region of interaction of this wave with the rear shock 
on the other side of the body, the shock strength has decayed suf- 
ficiently that the vorticity is of order 



(A. 10) 


for the planar case, whereas it is of order 

K\y 

o 


(A. 11) 
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only in the third approximation in the corresponding rotational flow. 
Thus, the velocity potential may be used to study this interaction, 
which appears as a second-order effect in this region. 

For flows about finite systems, the Lighthill region is even 
more distant, with a correspondingly greater decrease in the vortic- 
ity. In this region it is of order 



(A. 12) 


whence it must be accounted for only in the ninth approximation. 
The vorticity in the second approximation to the corresponding 
rotational flow is of order 


(A. 13) 
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Appendix B 


An Example : The Ci rcular-arc Airfoil 


We consider the solution for the wave structure emanating from 
a symmetric, biconvex airfoil/ of circular arc section. The body 
shape function is then 


£( ho - 0 i 
_ l-4tf 
- z7 > 

= O 


- >1 <. g fe <. , (B.01) 

^ < , 


and only the flow above the airfoil will be considered, as the flow 
is symmetric about the axis ^ O . 

From equation (3.35) we have that =. O in the regions de- 
fined by the vertical characteristics for / > ^ > and 


= 2.(f + 4vj) Efc - ^2. (B. 02) 

in the region swept out by characteristics from the body where, 
from (3.34) 


- 


g 

I ■+ 4v^ 


) 


whence 



3*1 y 

2 - 


(B. 03) 


(B. 04) 


These characteristics form a fan- like pattern and overlap the ver- 
tical characteristics from ahead and behind the body. Shocks must 
be inserted to render the solutions unique, and, by (3.12) must obey 


C/Sp _ 2 g 


(B. 05) 
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whence 



(B. 06) 


where the upper and lower signs refer to the leading and trailing 
shocks respectively. It is seen that the first-order wave pattern 
is symmetric about the vertical axis ( T = O ) for this symmetric 
(fore and aft) profile. The strengths of the shocks are given as 
a function of ^ by 


£6 r 1 = 



(B . 07) 


and are seen to vary as the inverse square root of the distance, 
for large distances, a well-known result for planar flows about 
general sections. 

The second-order solution is also zero along the freestream 
characteristics upstream of the leading shock, and is by (3.42) 
and (B.03) 


0, “ 64 (" ^ + 

+ + j 
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_ -J, 
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(B. 08) 


in the region between the shocks, and 


/ _ /•_£ [M z . -i— _ / f*! 1 2V\ 
^ ( 5 3/S*- + 3 )~j (+4^ + 3 / 


(B . 09) 


behind the trailing shock. The second-order shock displacement is, 
by (3.18) 




, HI 1, l)\ . — 3- _L TN 1 * 

- ( J2p + 6 Jv7+4^ + / + 4/5* 
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(B.10) 


for either the leading or trailing shock. The second-order incre- 
ment to the strength of the shock is thus 
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(B. 11) 


+ (^-"p + ^ ■ + "V) r4vj + (- ^ 


which, for large distances behaves as 



(B - 12) 


and is seen to increase the strength of the leading shock, and de- 
crease that of the trailing shock (since the coefficient in paren- 
theses is positive for typical gases) . 

Behind the trailing shock, the potential is a function of 
alone, whence the pressure in this region is zero to second order. 
That is not zero, but a function of , hints that there 

is a wave of the right-running (or returning) family that is locally 
of third order in strength, but has a cumulative effect of second 
order. The behavior of this wave cannot be calculated using the 
velocity potential, however, since the velocity and entropy pertur- 
bations are of the same order in this region. 

This wave is of interest because it must be included to give 
results consistent to second order for certain global integrals 
relating to lift and drag of the airfoil, and also because it even- 
tually intersects the rear shock on the other side of the airfoil, 
altering its position to second order. (In the symmetric case under 
consideration, the wave may be imagined to be reflected at the line 
of symmetry (the Jf -axis) and thus eventually interact with the 
upper rear shock.) See figure B.01. 

This interaction occurs only at very great distances from the 
airfoil, and at such distances that the shock has decayed suffi- 
ciently for a velocity potential to again exist to the order re- 
quired to study the interaction. Thus, a potential analysis may be 
used to study the interaction, but vorticity must be taken into 
account in determining the initial strength of the wave. 

The characteristic equations for steady, planar, inviscid (but 
rotational) supersonic flow are that 

along 


(B. 13) 


(B. 14) 
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Figure B.01: Wave Systems associated with Body in 

Supersonic Flight. 


Substituting expansions in terms of our wave structure variables 
find in the region behind the trailing shock, that the quantity 



/ 3 ^7 \ 7 *^ 


(B. 15) 


is constant along the lines 



The invariant (B.15) may be calculated at the rear shock and is, 
for our example. 



(B. 17) 


which, upon reflection from the J'-axis, gives a wave of strength 


0W 4 . n? 2-0 1 L „ 
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(B. 18) 


which is propagated along lines of constant *5 (linearized charac- 
teristics) . 

Now, the region in which this wave interacts with the rear 
shock is so distant from the airfoil that the shock has decayed to 
the point where a potential again exists; i.e., the entropy per- 
turbations are again small compared to the velocity perturbations 
in the wave (B.18). (See Appendix A) In this region (which we 
shall call the Lighthill region, after the discoverer of its im- 
portance (1954)), we define the new independent variables 


whence 




(B. 19) 



(B . 20) 


where the jump in potential is to be calculated from (B.09) and 
(B.18). From the latter we have 



(B. 21) 


or 


sp** 2 s)jJ7Zf “ ( 3 V* “aO , CB. 22 ) 

where the constant has been determined by matching the shock shape 
in this region with that in the wave structure region, using the 
asymptotic principle. This procedure results in 
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(B. 23) 




f JL _rn l £}) ?==- 

(sM z 3p + 3 


for the second-order correction to the 
shock. For very, very large distances 
becomes 


location of the trailing 
(i.e., large ), this 



(B. 24) 


This last result differs from the expression given by Light- 
hill (1954) for the analogous quantity in his analysis (equation 
(4.49)). The quantity given there seems to correspond to the shift 
in the rear shock beyond that calculated earlier by Friedrichs, and 
thus must be added to the original coefficient K (defined in his 
text between equations (4.09) and (4.10)). Also, since only the rear 
shock is affected (a fact noted earlier in the same article by 
Li gh thill) this quantity must be multiplied by the ratio of total 
mass flux excess in the simple wave to that in the rear portion 
(coming from downstream of the peak value of the airfoil). Thus, 
if for the upper shock LighthilPs equation (4.49) is replaced by 



our results are brought into agreement. 


Appendix C 


Mass^ and Momentum Fluxes in the Wave Systems 

That the theory presented here gives a consistent representation 
of the flow at large distances from the body may be verified by con- 
sidering the asymptotic form of the solution as it relates to the total 
forces on the body and to conservation of certain global integrals 
(such as mass flux).. Since the broad third-order wave behind the body 
(in the Lighthill region) has a cumulative effect of second order, this 
must be included in the analysis. We shall again restrict ourselves to 
planar flows, where the complete second-order solution may be obtained, 
and the strength of the third-order wave behind the trailing shocks may 
be determined. 

We shall consider integrating mass and momentum fluxes around con- 
tours sufficiently far removed from the body that the asymptotic forms 
of the solutions may be used, and show how these integrals correspond 
to conservation of mass and momentum, thus giving the lift and drag 
forces on the body. Consider the control volume shown in figure C.01, 
around which we will integrate mass and momentum fluxes. 



There will be no contribution to the flux integrals from a-b or e-f 
since the flow is parallel to the control surfaces; there will be 
a contribution from the difference of fluxes across c-d and f-a due 
to the increased entropy behind the wave systems; and there will be 
a contribution from the integration across the primary wave systems 
b-c and d-e. The effect of the third-order wave behind the body 
will be accounted for as it alters the position of the trailing 
shocks, and hence the contribution from the primary wave systems. 

For the symmetric circular-arc airfoil, we need consider only 
the upper half-plane. The excess mass flux in the primary wave 
system is 

(C.01) 



Thus, since 


and 




we have 


^ (UUcMt 1 / _± 4V\ 


(C . 02) 


(C . 03) 


(C.04) 


The excess mass flux in the entropy wake region, c-d, over that in 
the freestream, f-a, is given by 



(c ' 05) 


where the surface has been taken far enough downstream that pressure 
perturbations are negligible. Since entropy is conserved along 
streamlines, the integral may be taken at the shocks, and using our 
expression for the entropy jump at a shock, (C.05) may be written 
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and is seen to exactly balance the excess mass flux in the primary 
wave system, thus confirming that mass is conserved in the solution. 
Note that in (C.03) the value of the potential in the Ligh thill 
region, (3.53), has been used. 

The drag of the airfoil must be equal to the resultant of the 
axial pressure force and the deficit of axial momentum flux through 
the control surface. In the upper primary wave system, this is 
given by # 


which, for our circular-arc airfoil is 

(C. 07) 


(C . 08) 

which approaches zero as \l+4y) for large ^ , 
tributes nothing to the drag. *The deficit of axial 
in the entropy wake is 

and thus con- 
momentum flux 

= /f u f «"£> j z S' <ij j 

(C. 09) 


where, again, the plane of integration is assumed sufficiently far 
downstream that pressure perturbations are negligible. Thus, 



which is exactly the drag on the upper surface of the airfoil which 
is obtained by integrating the surface pressure times the axial 
component of the normal to the surface. 

The airfoil treated in the preceeding paragraphs has no net 
lift, since it is symmetric. To investigate lift forces, we con- 
sider an airfoil with the same circular-arc upper surface, but 
with a flat lower surface. Thus there will be no wave system from 
the lower surface, and the airfoil will have a net (negative) second- 
order lift. The control surface for this airfoil is shown in figure 
C.02; note that we have included integration across the third-order 
tail wave, which now is propagated to infinity below the section. 
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Figure C.02: Control Volume for Lifting Airfoil 



Here, the lift force must be balanced exactly by the resultant 
of vertical pressure forces and momentum fluxes, integrated around 
the control volume. We need consider only the integration across 
the primary wave system, b-c, and the reflected (Lighthill) wave 
system, d-e, since elsewhere the contributions are negligible. The 
contribution due to integration across the primary wave system is 



which for our case is 

* - 




?} 


(C.12) 
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In this case, there is no "Lighthill region" effect upon the position 
of the trailing shock, since there is no reflected wave coming from 
below the airfoil. Integration across the third-order wave coming 
from the upper, primary system 



which for our case is 


(C. 13) 



* \3 M 1 - 2(3 X + 3 f • 


(C.14) 


The sum of these two integrated pressure forces is thus 



(C. 15) 


which is the familiar Busemann result for the second-order lift of 
the section considered. 
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Appendix D 


The Structure of Shock Origin 


We consider the second-order structure near the point of forma- 
tion of a shock within the flow field, and its effect upon the point 
at which the shock is first formed. Since the characteristics are 
not revised after the first approximation, we cannot hope for the 
second-order theory to give detailed local information at the point 
of formation; however, we would hope that the solution at least ap- 
proximates any such local behavior in some global sense. Such is 
the case, and we will, in fact, show that the shock displacement pre- 
dicted by the second-order theory at the point of formation corres- 
ponds exactly to the shift in virtual origin the shock must have to 
account for this local behavior. 

To demonstrate this fact, we consider the planar flow over a 
parabolic compression-expansion bump. The first-order theory pre- 
dicts formation of a shock of finite initial strength for this case 
within the flow field. (See Figure D.01) 



Figure D.01: Geometry of Shock Origin Problem. 
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First, the wave structure solution will be found, and the second- 
order displacement of the shocl^ at its point of formation will be 
calculated. The geometry of the actual characteristics (accurate 
to second order) in the vicinity of the point of formation will 
then be considered to show how this displacement corresponds to 
the true local behavior. 

1 . S econd-o rder wave structure solution 

We consider the planar flow past the curved wall described by 


■P(r b -) = ° , ‘S'b < o 

- r b z , o < < l ) 

= 2.- ( 2 -g b f > i < g b < e , 
= 2. , a < ■eu . 

Then, we have by (3.35) that for 0<t b < 1 , 



-g* 
1 - 2 . ‘ 


and for 1 < Z , 


& - - 2 ♦ 


(z-wt 

1 + 2*1 


(D .01) 


(D . 02) 


(D . 03) 


A shock starts at (0,^.) with finite strength and obeys 


£<>= . 


(D . 04) 


The second-order solution, by equation (.3.42) is seen to be, 
for 0<- g b < 7 (and ) 


. 4TM\ 3A ) .. » ^ -2- _ grAf*) — ; 

$ = ( 3 3jS*- + 3 4 ( 3 + M* 3/S* 

and for I < < 2. , 


(D. 05) 


$ = f-?r- 3^*+ C 2 - K) 2 * (§- " 4) 

m ^ 9 > * /T~\ 


+ (-3 ■ ^(^-0(2- SJ 

, 44TM x ; 

+ 3 3/S*- 


(D . 06) 
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Evaluating this along the position of the first-order shock yields 





CD. 07) 


and since 



CD. 08) . 


the second-order shock displacement is given by (3.18) as 



(D . 09) 


Finally, the displacement due to second-order effects at the 
point of formation is 


(*) - 


2 .-+ 


j~ zn$ l 

3f&*- 



(D.10) 


2. Interpretation in terms of local structure 


The effect of second-order corrections to the characteristics 
is to give a local structure to the point at which all the charac- 
teristics from the compressive portion of the body converged accord- 
ing to the first-order theory. The geometry of these characteris- 
tics (accurate to second-order) and the shock wave in the neighbor- 
hood of the point of shock formation is now considered. The left- 
running characteristics are given in terms of our wave structure 
perturbation potential as 




CD. 11) 
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Since the flow may be described as a simple wave to second 
order, the flow deflection angle 

e= -T(f> w + -p + ••• , (D.12) 

and, hence, also 0^ , are constant along these characteristics. 
Application of the boundary condition at the body surface gives 

4 = - fk) * t(- M f f “) - §$("'} , (D. 13) 


whence the equation of the characteristics becomes 


•. tp.uj 


Thus, for the parabolic bend described by 


£(e 0 = sj- , 


CD. 153 


we have 




CD. 16) 


which may be integrated to give the equation of the family of char- 
acteristics 


°> CD,17) 

accurate to second order. Detailed information about the charac- 
teristics from only the compressive portion of the body is essen- 
tial in describing the early stages of formation of the shock, as 
will later be seen. 

The shock will first be formed in the vicinity of the envelope 
of these characteristics. This envelope is determined by elimin- 
ating between the preceeding equation and the equation obtained 
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by differentiating it with respect to the parameter, : 

-1+2^- = O. CD.183 

Thus, the envelope is parametrically described by . 


£ = 


1 


- 4 2 a) *£ , 


(D.19) 


from which it is seen to be the upper half of the parabola with 
vertex at ( 0^2-) and opening to the left if the factor 


> 0 _ r . JL zrM* 9 

^ ^ M 7 - {&■ * 5 ^ 


(D. 20) 


is positive, or the reflection about the origin (vertex) of that 
curve if K is negative. In the vicinity of the origin of the 
shock, then, the problem naturally divides itself into two cases 
depending upon the algebraic sign of fc. . (For air, treated as 
a calorically perfect gas, O when M = 1.33, is negative for 
lower supersonic Mach numbers, and is positive for greater Mach 
numbers.) For positive X , the shock starts with zero strength 
at .( 0 > Yz> ) -- the same point at which the finite-strength, first- 
order theory shock originated -- and builds up to full strength in 
a region of scale 0 (t) in the plane. For negative ^ , 

the shock again starts with zero strength at the point on the en- 
velope corresponding to the characteristic from the steepest por- 
tion of the body, and develops for a short distance (again of 
0 (t) in the plane) with disturbed flow on both sides. See 
figures D.02 and D.03. 

If ^=0, the envelope collapses to the point ( 0> ) as 

in the first-order theory, and the shock starts from that point 
with finite strength. Since the total length scale in this local 
development is short, in considering the path of the shock, we 
will use simply the characteristic -bisector rule with negligible 
error. 
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a. 


Case I, *> 0 : 


For the case when X is positive, the shock starts at ( Oy^Z-) 
with zero strength, and its path is determined as the bisector of 
the freestream characteristics and those coming from the body for 
O <■ Eb <. 1 . The shock will intercept these characteristics 

coming from the body before they can form the envelope calculated 
in the preceding section. (See Figure D.02). We set up a local 
coordinate system with origin at the point ( O, Vz. ) and 



(D. 21) 


by 


Then, the characteristics coming from the body are described 


* *-2E 


SfeJ > 0 < £ lo < 1 > (D . 22) 


and the slope of the shock at every point is determined by the bi- 
sector rule 



(D. 23) 


Eliminating from this last equation using (D.22) gives 


<hL _ 



(D • 24) 


This equation cannot easily be solved in general due to its non- 
linear nature, but the unique solution passing through the origin 
is easily verified to be 


X =- 


_3_ 

4 - 


3 


(D . 25) 


The shock follows this curve until the point where it intersects 
the characteristic coming from Ew-1, (i.e., from the inflection 
point on the body) . This point is found to be 
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X = - *3 




<j ■= ^3 • CD. 26) 


From this point on, the shock is determined to satisfy the appropri- 
ate relations between the freestream characteristics and those com- 
ing from the expansive portion of the body, as in the wave structure 
solution. If we analytically continue this shock shape back to the 
value yf st y'z- Ci.e.» y = 0) , which is the origin according to first- 
order theory, we find 


*/ 



or equivalently. 


¥ > ( D - 27 ) 

which is precisely the displacement predicted by the second-order 
wave structure theory. Thus, the displacement of the origin of 
the shock predicted by the second-order wave structure theory cor- 
responds to the displacement of the virtual origin the shock must 
have to account for the local formation process. 

b . Case II , tt<- O : 

For the case when is negative, the shock starts at the 
point on the envelope of characteristics corresponding to i 

with zero strength. The path of the shock is then determined as 
the bisector of the characteristics coming from the compressive 
( O < t <. 1 ) and expansive ( 1 < < 2L ) portions of 

the body, respectively. (See Figure D.03). The shock will again 
intercept the characteristics coming from the forward portion of 
the body before the envelope is formed. Using the same local co- 
ordinate system as in Case I, the shock slope obeys 

% * tD - 28) 

when the body has been assumed straight for with no loss in 

accuracy. Again, a general solution to this equation is not found, 
but the particular solution satisfying the appropriate initial con- 
ditions is 






% 


CD. 29) 
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The point after which the shock again becomes determined by the 
freestream characteristics and those coming from the expansive por- 
tion of the body is 


X =* 

u = - ha 


(D. 30) 


For this case, since the shock is formed for v j<^z. , we merely 
the value of Jg^ when i ^z. • It is ' 


- _ _ 2TM\ 2m 

^ + M l 3S 7 - + 3 > 


(D.31) 


again corresponding to the displacement of the virtual origin of 
the shock predicted by the second-order wave structure theory. 
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Appendix E 


Some N ume rical Comparisons 


The primary importance of higher-order asymptotic theories 
is the justification they provide for the lower-order theories, 
by showing that they are the first step(s) in a rational approx- 
imation to the solution in some limit. Since the series thus 
produced are not necessarily convergent, the inclusion of higher- 
order terms can worsen numerical accuracy for calculations at 
finite values of the appropriate small parameter -- a fact for 
which asymptotic series are somewhat notorious. It is, however, 
still a matter of interest to see how the inclusion of such higher- 
order terms affects numerical results, and, for that reason two 
such examples are included here. The two simplest examples in 
each, planar and axisymmetric flows, are considered because (1) 
exact inviscid calculations were readily available for comparison, 
and (2) only the simplest geometries are conducive to solution in 
the axisymmetric case. The shock angles predicted for these flows 
are calculated and compared with the exact inviscid solutions in 
the following sections. 

E.l. Semi-infinite Wedge (Planar flow) 

The body shape function for a semi- infinite wedge of semi- 
angle -r is 


•f(Eb) - O > 


*b< O , 

O < ‘Et • 


(E. 01) 


Thus, in the region of disturbed flow, the first-order potential 
is 


^o r = - / > (E . 02) 

and the shock separating this region from the undisturbed stream 
obeys 



(E . 03) 


The second-order potential, evaluated at the first-order 
shock location, , is then 
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(E.04) 


^Cc oi “ ( f + 2M* ~ ^ ~ > 

whence 


/ . _!_ 3TM 7 - Us -J, \ 

— ( 1 + 2M 2 - " Sp* + I'Z “ / iz) v ^ (E . 05) 

The shock location, to second order is thus 

** = {l'*p( Z +M z - ^”4)J , CE - 06) 

whence 

coiS sl = (e. 07) 

Numerical values for the shock angle for flows in air (treated 
as a calorically perfect gas with P = 1.2, ^ = -1.1, and 

V =0.4) are compared with results of exact inviscid calculations 
from the Ames Research Staff (1953) in Figures E.01 and E.02. The 
second-order theory is seen to give considerable improvement in 
numerical accuracy even for quite high Mach numbers (e.g. , = 5). 

E.2. Semi-infinite Cone (Axisymmetric flow) 

The calculation for the cone is much more complicated than the 
wedge for two reasons. First, an additional scaling region near the 
body must be considered, and this "slender body" solution used to 
determine the local solution by matching. Secondly, both of these 
inner solutions must be determined to third order to determine the 
true second-order wave structure solution. 

All of the details of this extremely tedious calculation will 
not be given here. The general nature of the slender body theory 
from the viewpoint of matched asymptotic expansions may be found 
in Ashley and Landahl (1965) . No peculiar difficulties are en- 
countered in extending the solution to third order (for the case of 
a cone!) beyond algebraic tedium. Only the radial derivative of 
the third-order solution was obtained, since this determines the 
local solution. (The complete solution is fiecessary if we are inter- 
ested in pressures near the body.) 

For the flow over a slender cone of semi-angle r , the first- 
order shock strength is of 0 (t 4 J, whence a velocity potential exists 
to the required third order (which describes perturbations of 0(tO 
in the local region. This local solution is of the form 
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mach number 


Figure E.02: Shock Angle on Sixteen-degree Wedge. 
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where the conical variable 




(E.09) 


has been introduced, and particular integral terms which do not 
enter directly into the determination of the wave structure so- 
lution have been omitted for clarity. Of these omitted terms, 
some will automatically match with particular integral terms in 
the wave structure solution, and the others will be negligible. 
The constants in this local solution are determined from the 
slender body matching to be 


» 

"tS = (2- 'M rjz ~ 1 ” " ~js*- > 

C + +^2-')} (E.10) 

+ - 1 - -tm 1 } 

v- rM v^. 

The large distance (i.e., 1 ) behavior of the local so- 

lution is given by 



which must match with the small-distance behavior of the wave 
structure solution, 



_ r - (? - ^ 5( «V i / h 

If + - 


(E . 12) 
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where, again terms which do not enter directly into the determina- 
tion of the unknown functions in the wave structure solution have 
been omitted. The wave structure solution is thus determined by 
the conditions 


Q(*u) 

Hteb) 


- 


(E. 13) 


where 


\/z% - y{t * } . 


(E.14) 


The first-order shock position is then found to be 



(E . 15) 


The one-and-one-half-order correction is 





(E . 16) 


And the second-order correction is 


r, = {-■£- ^4^ 4 ("z 1 - 3 

+2.- izo 1 ^- % - ilftVe?} 1 ] 1 - 

Thus, the shock angle is given by 

cot « /3^ 1 - i 5^{t + 2MVB + 

+ ((0-3 + 


(E. 17) 


CE. 18) 


Numerical values for the shock angle for flows in air are 
again compared with those from exact inviscid calculations (Ames 
Research Staff (1953)), and are presented in Figures E.03 and 
E. 04. 
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Figure E.04: Shock Angle on Fifteen-degree Cone. 


74 


Note that, generally, the one-and-one-half-order solution does not 
give results that are appreciably better than the first-order theory, 
and we must go to the full second-order theory to achieve noticeable 
improvement. That the second-order theory, which is based upon the 
local third-order theory gives increased accuracy for these finite 
values of body thickness, is probably due in part to the fact that 
the theory is of second order near the shock. Third-order theories 
are often very bad for numerical calculations at finite values of 
the asymptotic parameter. 
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